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Abstract 

We show that the reported anomalies in Rf, and Rc can be interpreted as the effect of a 
heavy vector boson V universally coupled to u- and d-type quarks separately and nearly de- 
coupled from leptons. This extra vector boson could then also naturally explain the apparent 
excess of the jet rate at large transverse momentum observed at CDF. 
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1 Introduction 



On the whole the electroweak (EW) precision tests performed at LEP, SLC and at the Tevatron 
have impressively confirmed the formidable accuracy of the Standard Model (SM) predictions. 
There are only a few hints of possible deviations and our hopes of finding new physics signals are 
confined to them. At LEP the observed values of Rb and Rc deviate from the SM predictions by 
about 3.5cr and 2.5a, respectively At CDF an excess of jets at large with respect to the QCD 
prediction has been reported 0. None of these observations provides a very compelling evidence 
for new physics as yet, because of the limited statistics and of possible residual experimental 
systematics. The Rf, value is relatively more established, in the sense that it was first announced 
in 1993 and is insofar supported by the analyses of all four LEP collaborations with several 
independent, in principle clean, tagging methods. From a speculative point of view it is not 
implausible to have a deviation in the third generation sector. Also, a moderate increase of Rb 
with respect to the SM (of roughly half of the present excess) would bring the value of a<j(M^) 
measured from the Z widths in even better agreement with lower energy determinations. The Rc 
evidence is much less believable both from the experimental and the theoretical points of view. In 
absolute terms it is a large deficit, that would overcompensate the Rb excess. Thus these results, 
if taken at face value, would demand a deviation from the SM in the light-quark widths as well, in 
order to reestablish the observed value of T^, which is measured with great experimental accuracy 
and agrees with the SM. After all, in this context charm is alike any other first or second generation 
quark, while beauty could be special, being connected to the heavy top. If one literally believes 
the data, then one must accept an accurate cancellation among the new physics contributions 
to light and heavy quarks. But the perfect agreement of the leptonic widths with the SM, up 
to a fraction of MeV, clearly poses the problem of how to naturally shift the light quark widths 
without affecting the leptonic ones as well. Finally, the significance of the CDF result on jets 
entirely depends on the calculation of the QCD predictions at large E^, which could to some 
extent be questioned. For example, it was recently pointed out @| that it is possible to shghtly 
increase the large-x gluon densities without deteriorating the standard overall fits to low energy 
data, and thus partly explain a large fraction of the high-ii^x jet discrepancy. 

All these words of caution being said, in this note we consider the challenging task of quan- 
titatively explaining in an admittedly ad hoc but relatively simple model all the three observed 
deviations discussed above. We introduce a heavy vector neutral resonance V, singlet with respect 
to the standard gauge group SU{3) x SU{2)l x U{1)y and with a mass in the TeV range. We 
allow this new resonance V to have a small mixing with the ordinary Z gauge boson, and therefore 
to contribute to the Z decays. We observe that while in the data 5{Rb + -Rc) is large and negative, 
5{?>Rb + 2i?c) is only about la away from zero. This suggests to take universal couplings of the 
V to the three generations of fermions separately for up, down and charged leptons. Since the 
leptonic width is in perfect agreement with the SM, the leptonic couplings of V must be much 
smaller than those needed for the quarks to explain the deviations observed via 5Rb and via 5Rc, 
and we shall take them as approximately vanishing (at a less phenomenological level, one must be 
prepared to add new, presumably very heavy, fermions to compensate the anomalies). Then the 
products of the amount of mixing (which is severely constrained by the data) times the couplings 
of the V to up- and down-type quarks are fixed by imposing that the observed values of 5Rb and 
of 5{?>Rb + 2Rc) be approximately reproduced. We have at our disposal five parameters to do that: 
the amount of mixing. My (that for a given mixing fixes ei = 5p), the left-handed coupling to 
the (m, d)i^ doublets, and the two right-handed couplings to the and c/r singlets. So the game 



would be trivial, were it not for the fact that couplings to quarks as large as those required by 5Rb 
and SRc would tend to produce too large effects in the distributions of large- i?T jets measured at 
the Tevatron. We can then adjust My > 1 TeV and the left and right couplings in such a way as 
to obtain a reasonable fit to both LEP and CDF anomalies, without violating, to our knowledge, 
any known experimental constraint. The details are given in what follows. 



2 Effects on LEP and SLC observables 

The tree level neutral current interaction can be written in terms of the unmixed interaction states 
Zq and Vo, coupled respectively to the ordinary standard model neutral current (J3L — sin^ OwJem) 
and to an additional current Jn- The vector and axial couplings of the gauge bosons Zq and Vq 
are defined by: 

z cos uw i 

The Zq couplings are the standard ones 

v's = nL-2sm^9wQ\ 0^3 = -Tk (2-2) 

where T^^ is the third component of the weak isospin of the fermion i, and its electric charge. 

We assume that the new gauge boson Vq couples only to the quarks and has zero (or negligible) 
couplings to the leptons. We also assume family-independent couplings. The new interactions can 
be then be expressed in terms of three parameters x, y„ and i/d'- 



(2.1) 



v"^ = x + yu 

d 



V 



N 



-X + Vd, 



(2.3) 



where the superscripts u and d refer to up-type and down-type quarks. 

In presence of a mixing, the mass eigenstates Z and V are given by a rotation of the unmixed 
states Zq and Vq: 



Z — cos^Zq -I- sin^V^) 
V — — sin^Zo-|-cos^Vb 



Due to the mixing, the p parameter, defined by 

P 



M|cos2 Ow 



(2.4) 



(2.5) 



receives a tree level contribution ApM, which in term of the V mass and the mixing angle ^ is 
given by: 



M 



- 1 



sin^ ^ ~ 



/Ml 



V 



(2.6) 



At LEPI the observables get corrections from the presence of V through the mixing with the 
ordinary Z and through the shift in the p parameter. Contributions from direct V exchange are 



negligible at the Z pole, but will be taken into account later on in our study of the Tevatron jet 
observables. 

The deviation of a LEPI observable, linearized in Ap^ and ^, can therefore be expressed as: 

— = AoApM + Boi. (2.7) 

The coefficients Aq are universal and depend only on the SM parameters and couplings, while 
Bq also depend on the Vq couplings f]y and a]y 0]. In Table I we give the numerical values of Aq 
and the expressions for Bq for the observables of interest, as functions of the parameters x, ?/„ and 
yd introduced in eq. (p73|). In Table I we also present the experimental data used in the present 
analysis fll], together with the Standard Model predictions p for rrito-p = 175 GeV, mu = 300 GeV 
and as{Mz) = 0.125. They include the one-loop electroweak radiative corrections. The Z mass 
was fixed at the experimental value Mz = 91.1887 GeV. 



Quantity 


A 


B 


Exp. values 


SM values 


Pull of the fit 


Tz 


1.36 


-0.92x - 0.49y„ + 0.37yd 


2496.3 ±3.2 


2497.4 


1.72 


Ri - Th/Ti 


0.34 


-1.31a; -0.70y„ + 0.52yd 


20.788 ±0.032 


20.782 


-0.35 


(Jh 


-0.030 


0.52a: + 0.28y„ - 0.21yd 


41.488 ±0.078 


41.451 


-0.32 


Rb = Tb/Fft 


-0.094 


-3.16x + 0.70y„ + 0.29yd 


0.2219 ±0.0017 


0.21569 


-1.41 


Rc = Tc/r/i 


0.12 


6.20a; - 1.43y„ - 0.59yd 


0.1543 ±0.0074 


0.17238 


1.62 


Mw/Mz 


0.71 





0.8802 ±0.0018 


0.8808 


0.94 


Ai 


18.50 





-0.15066 ±0.00276 


-0.14334 


0.98 


Ab 


0.23 


-0.31a;- 1.72yd 


-0.841 ±0.053 


-0.9342 


-1.79 




1.70 


2.37a; + 5.36yu 


-0.606 ±0.090 


-0.6662 


-1.00 


^FB 


18.15 


-0.31a;- 1.72yd 


0.0999 ±0.0031 


0.10042 


1.20 


^FB 


19.63 


2.37a; + 5.36y„ 


0.0725 ±0.0058 


0.07161 


0.76 



Table I : Coefficients A and B, defined in eq. (|2.7|) , for various electroweak observables, 
together with their experimental values and SM theoretical predictions for rritop = 175 GeV, 
rriH = 300 GeV and as{Mz) = 0.125. The corresponding is equal to 26.73. In the last column 
we report the pull values ((fit-exp)/(j) for the final V fit with x = —1, ?/„ = 2.2, yd = and 
^ = 3.8 ■ 10~^. The iii this case equals 14.72. 

The deviations in Table I are computed from the tree level formulas for the partial widths 



and for the asymmetries 



G 

nZ - //) = -^PNc [{vi^ff + [ai^ff] , (2.8) 



The forward-backward asymmetries are given by: 



(2.10) 



In eq. (W^ ) Nc = 3 for quarks and Nc = I for leptons, and in eq. (|2.8| ) and (p.9|) the effective 
vector and axial- vector coupling v^jj and a{jj are superpositions of the corresponding Zq and Vq 
couplings: 



a 



f 

eff 



COS C, ag + sin ^ a{f (2-11) 



In computing the deviations due to the new vector resonance V, it is sufficient to consider the 
tree level expressions for the observables, because the corrections are proportional to ApM or ^, 
that are both constrained to be quite small (of the order 10~^) by the current electroweak data. 

We keep fixed the input parameters a, Gp, Mz, and take into account the modification of the 
effective Weinberg angle given in eq. because of the shift in the p parameter. One finds 

d{sm Ow) = ^P^i (2-12) 

cos"' 20w 

The loop effects due to the heavy gauge boson V are quite small and we will neglect them. 



3 Fit to the LEP and SLC data 

In this Section we constrain the free parameters of our extended gauge model by performing a fit 
of the eleven independent observables of Table I. The parameter space of the model includes the 
couplings X, i/u, Vd and the two parameters ^ and My. ApM is related to the previous parameters 
by eq. (H). 

We have minimized the function keeping My fixed at different values: it turns out that the 
best fit central value for ApM stays almost fixed, by varying My, at the value 

ApM ^ 0.0011 (3.13) 

This implies, from eqs. ( p.6|) , that the mixing angle ^ decreases with My: 



^ ~ VO.OOll^ (3.14) 

The parameters x, yu, and yd are multiplied by the mixing angle C, in the expression ( p.7| ) for 
the deviations: that means, from eq. (|3.14|) , that their best fit values will scale with My as 1/^, 
i.e. 

Mz 

The fit, for a choice My = 1000 GeV, leaving the four parameters ^, x, yu, and yd as free, 
gives: 

e = (3.0^?:^) ■ 10-^ , x = -iAt\-l 

yu = 5.3+t-l , Z/d = 2.911:6 (3-16) 

We have quoted the standard errors, corresponding to = Xmm + 1- The fit is weakly sensitive 
to the parameter yd, as the large error indicates. We take advantage of this by constraining the 
fit with yd = 0. The other parameters of the fit with yd = turn out to be: 

e = (2.8in) ■ 10-' , X = -2.11:11 , yu = ^.5ttt (3.17) 



Figure 1: 70% confidence level ellipsis in the plane yu versus x, for My = 1000 GeV, 
yd = and ^ = 2.8 • 10^^ (left figure) or ^ = 3.8 • 10"'^ (right figure). 



where we have fixed, as before, My = 1000 GeV. Here = 11.41. For other values of My, 
the scahng formulas eq. ( |3.14 ) and ( p. 15 ) are an excellent approximation. The central values 
correspond to quite large couplings that would be incompatible with the CDF data, as shown in 
the next Section. 

The parameters in eq. ( p.l7| ) are strongly correlated. The correlation between the parameters x 
and i/u is easily understood once noticed, from Table I, that the ratio of the coefficients multiplying 
X and i/u in the formulas for the deviations is the same, 1.87, in the observables Tz, Ri and ah- The 
relative high precision data are in excellent agreement with the SM predictions, and this induces 
a strong anticorrelation between the two parameters. 

In fig. |l|a we plot the 70% confidence level ellipsis in the plane ?/„ versus x, keeping ^ fixed at 
the best-fit value 2.8 ■ 10~^. From the figure, one can see that at this confidence level the points 
closest to the origin are at the position x ~ —1.3, yu — 2.8. 

In fig. we present the analogous ellipsis for the higher value ^ = 3.8 • 10^'^. Increasing ^, the 
elliptical region moves toward the origin, because the higher mixing angle forces the parameters 
X, Uu to smaller values. For the value of ^ = 3.8 ■ 10~^ the closest points to the origin are located 
at X ~ —1.0, Uu — 2.2. Moving away from the best-fit value of ^, the value increases: for 
^ = 3.8 ■ 10~^, X = —1.0, and yu = 2.2, one obtains = xLin + 3-3, still in the 70% confidence 
level region of the three parameters fit of eq. ( ^.17] ). 

In the last column of Table I we quote the pull values, given by (fit - exp)/cr, for ^ = 3.8 ■ 10~^, 
X = —1, Uu = 2.2 and yd = 0: the discrepancies in Rb and Rc are reduced. We stress again that 
these are not the best fit values for the parameters, but they lead to an effect on jet observables 
which is quite compatible with the CDF observations, as we shall now discuss. 



4 Comparison with the Tevatron jet distributions 

A vector resonance V with such large couphngs as obtained from the fits of the previous Section 
is hable to produce visible effects in hadronic collisions. There it can be directly produced via the 
Drell-Yan mechanism if the mass is not too large, or can lead to effective interactions between 
quarks via virtual exchange. The net result is a growth of the inclusive Et distribution of jets at 
large E^, relative to the standard QCD expectations. Using the V^couplings defined in eq. (|2.3|) , 
it is easy to evaluate the following quark-quark scattering amplitudes (amplitudes for crossed 
channels can be easily obtained from these ones): 

1 — 1 — 



(47r)2 {Att] 

-a,a s n{ ^^^^ _ ^2 ) ^ ^^^^^^ + ^^^^ _ ^2 ) + ^MvTv' ^^"^ ^ ^ 



'-IV/ J- v\ 

1 



2 1 1 

3 (t - M^) + iMyTv {u - M^) - iMvTv^ ^ 

2 2 

2^'^?[(t_M2)2 + M2r/ + («-M2)2 + M2r/^-' ^^-^^^ 

1 Y^\A{qq' qq')\^ = J—Y^\AQ^^{qq' qq')\' + 



(47r)2 ^' (47r) 



where Aqc^i is the standard QCD amplitude, IZ denotes the real part, a' = g"^ / {IGn cos"^ 9w) ~ 
0.010 and Ty is the V total decay width given by: 

Tv = 2Nga'Mv{2x' + yl + yj) (4.20) 

Taking the My — > oo limit one recovers^ the standard results obtained in presence of an effective 
4-quark coupling . 

Fig. H shows the deviations induced by the couplings to the V on the jet inclusive E^ distri- 
bution at the Tevatron. The quantity: 

is plotted as a function of jet E^ for different values of x and chosen in the range favoured by 
the EW fits. This is compared to the CDF data 0, represented in the figure as: 

^"'"'/"^^ - 1 (4.22) 

The calculation of the V contribution incorporates the full set of QCD processes, including re- 
actions initiated by gg and gq. Only LO diagrams are considered, as no NLO calculation for 

■^Up to some misprint contained in the standard literature. 




Figure 2: The effect of V exchange on the inclusive distribution of jets at the 
Tevatron. The different curves correspond to increasing values ofuu, from 2 to 4- The four 
displays correspond to x = —0.5, -1, -1.5 and -2. The quantity {C DF — QC D) / QC D ^ 
is shown by the points. We used dashed or continuous lines to indicate whether the V 
width is smaller or larger than 500 GeV. 




Figure 3: The effect of V exchange on the invariant mass distribution of di-jet events 
at the Tevatron. The different curves correspond to increasing values of y„ from 2 to 
4- The four displays correspond to x = —0.5, -1, -1.5 and -2. The quantity {CDF — 
QCD)/QCD ^ is shown by the points. 



the l^-exchange contribution is available. The calculation was performed using the MRSA set of 
parton densities, and a renormalization scale /i = E't- We verified that the quantity displayed 
in fig. ^ is very stable under changes of these parameters. We also expect that NLO corrections 
should not affect significantly our results. 

As the figure shows, the extreme choice x = —1, = 2.5 allowed by the EW fits is fully 
consistent with the CDF data. A similar conclusion can be reached by examining the di-jet mass 
distribution, shown in fig. ^. Notice that the peak structure disappears for too large couplings, as 
the convolution of the large width and the falling parton luminosities smears away the resonance. 

5 Conclusions 

Deviations from the SM in Rh, R^, and in CDF jets have been reported. They do not yet constitute 
compelling evidences for new physics. Nevertheless one may want to take them at their face values 
and look for some new effect to explain them. We introduce, as a simplest object, a new heavy 
singlet vector boson, with some mixing to the Z and direct couplings to quarks, the same for all 
up and the same for all down quarks, we perform the overall fit to LEP data, and see whether we 
can also explain CDF jets. This is possible , within the errors, with a vector boson of mass larger 
or of the order of 1 TeV, weakly mixed to the Z, but rather strongly coupled to the quarks. We 
do not attempt at this stage any deeper theoretical construction. 

After completing this work we received a paper where similar ideas are discussed 

6 ADDENDUM: Low energy neutral-current data 

The data analyzed in the main body of this work do not include low-energy neutral current 
experiments. The present Addendum is devoted to size the impact of deep inelastic neutrino 
scattering on the allowed region in the parameter space. 

The relevant information is contained in table II, where, with the same notations used above, 
we list experimental data, SM expectations and deviations for the four parameters gj^ ^ and Ol^fi 
characterizing z/-hadron scattering 



Quantity 


A 


B 


Exp. values 1^ 


SM values 


Pull of the fit 


9l 
9\ 
Ol 

0R 


2.71 
-0.60 
-0.07 
0.0 


-0.45 X 
-9.33yu + 4.67yrf 

-1.04a; 
0.50?;„ + 1.00yd 


0.3017 ±0.0033 
0.0326 ±0.0033 
2.50 ±0.035 

4 fiS+O-^^ 

^■^°-0.28 


0.303 
0.030 
2.46 
5.18 


0.76 
-1.63 
-0.79 
1.64 



Table II : Coefficients A and B, defined as in eq. (2.7), for low-energy neutral current 
observables, together with their experimental values and SM theoretical predictions for rritop = 
175 GeV, rriH = 300 GeV. In the last column we report the pull values ((fit-exp)/cr) for x = —2.1, 
Vu = 4.3, = and e = 2.3 ■ lO'^. 

Including also the four low energy observables in the fit, fixing as before My = 1000 GeV 
and yd = (the fit does not improve significantly releasing this parameter) and leaving the three 



parameters ^, x and i/u free to vary, one obtains: 



e = {2.3t\i) ■ 10 



X 



-2 1+^-° 

'^■'--3.7 5 



Vu 



4.3 ±2.9 



(6.23) 



The of the fit is 20.2, while the SM, for the values listed in table II, gives = 30.7. We recall 
that, by omitting the low-energy data, we obtained: 



e = (2i 



h0.9N 
-I.4J 



■ 10" 



X 



9 1 +0.8 

-Z.i_2,4 , 



Vu 



4 ^+^-'^ 



(6.24) 



Comparing (|6.23|) with ( |6.24| ), one notices that the low-energy data do not affect the results of the 
fit in any significant way: central values and errors are essentially determined by the LEP data 
alone. 

As we have discussed in the main body of the work, the central values in ( |6.23D or ( |6.24D give a 
too strong enhancement in the inclusive jet cross section at large Et, incompatible with the CDF 
data. In the low region, the values ^ ~ 3.8 ■ 10"'^, x ~ —1.0 and i/u — 2.2 previously retained 
remain a good compromise also when including in the fit the set of low energy data, which as we 
have shown do not practically influence our analysis. 

We have also included, in a following step, the weak charge Qw of Cesium |]10| measured in 
atomic parity violation experiments: the result is a small (~ 10%) decrease of the central values 
of the parameters x and i/u in (|3.17| ). 



In conclusion the situation remains practically unchanged after inclusion in the fit of low energy 
data and we hope that future high energy data will clarify the problem. 



Acknowledgements 

We thank Alain Blondel for the suggestion to include the neutrino deep inelastic data in our 
analysis, Paul Langacker for useful comments and Kevin McFarland for stimulating criticisms. 

References 

[1] The LEP Collaborations Aleph, Delphi, L3, Opal and the LEP Electroweak Working Group, 
CERN-PPE/95-172. 

[2] F. Abe et al., CDF Coll., FNAL-Pub-96/020-E. 

[3] J. Huston et al., Michigan State Preprint MSU-HEP-50812, |hep-ph/9511386i 



[4] G. Altarelli, R.Casalbuoni, S. De Curtis, F.Feruglio and R.Gatto, Mod. Phys. Lett. A 5 
(1990) 495; Nucl. Phys. B 342 (1990) 15. 

[5] D. Bardin et al., in Reports of the Working Group on Precision Calculations for the Z Reso- 
nance, eds. D. Bardin, W. Hollik and G. Passarino, CERN Yellow Book 95-03, p. 7. 

[6] M.A. Abolins et al.. In Snowmass 1982, Proceedings, Elementary Particle Physics and Future 
Facilities, 274-287; 

E. Eichten, K. Lane and M. Peskin, Phys.Rev.Lett. 50 (1983) 811. 



[7] E. Buckley-Geer, for the CDF Collaboration, FERMILAB-CONF-95-316-E, Sep 1995. 4pp. 
Presented at International Europhysics Conference on High Energy Physics (HEP 95), Brus- 
sels, Belgium, 27 Jul - 2 Aug 1995. 

[8] P. Chiappetta, J. Layssac, F.M. Renard and C. Verzegnassi, PM/96-05, piep-ph/960130€ 

[9] L. Montanet et al. Review of Particle Properties, Phys. Rev. D50(1994) 1173. 

[10] M.C. Noecker, B.P.Masterson and C.E. Wieman, Phys. Rev. Lett. 61(1988) 310; V.A. Dzuba, 
V.V.Flambaum and O.P. Sushkov, Phys. Lett. A 141(1989) 147 ; S.A. Blundell,W.R. Johnson 
and J.Sapirstein, Phys. Rev. Lett. 65 (1990) 1411. 



